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Abstract 

One loop corrections to the muon decay are studied in a popular and self-consistent 
version of the Left-Right symmetric model. It is shown quantitatively, that the cor- 
rections do not split into those that come from the Standard Model sector, and 
some decoupling terms. For a heavy Spontaneous Symmetry Breaking (SSB) scale 
of the order of a least 1 TeV, the contributions from the top quark have a loga- 
rithmic behaviour and there is a strong quadratic dependence on the heavy Higgs 
scalar masses. The dependence on the light Higgs boson mass is small. The heavy 
neutrinos are shown to play an important role, although secondary in comparison 
with the heavy scalar particles as long as the heavy neutrinos' Majorana Yukawa 
coupling matrix Hm obeys unitarity bounds. 



1 Introduction 

Embedding tlie Standard Model (SM) into a larger gauge group increases the 
number of degrees of freedom. For the Left-Right Symmetric Models (LRSM) 
based on the SU (2) SU (2) r®U (1) b-l gauge group [1,2] these are connected 
with new fields and interactions. The model is complex with extra particles 
of different types. New neutral leptons, charged and neutral gauge bosons, 
neutral and charged Higgs scalars appear. There are many different versions 
of the LR models with equal or different left and right gauge couplings gL,R, 
and specific Higgs sector representations. This robust structure is a chalenge 
and a good theoretical laboratory for testing many phenomena beyond the 
SM. 
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The purpose of the present work is to study numerically one loop corrections 
to muon decay which come from the extended gauge sector of the LRSM. Since 
the history of the LRSM is already quite long, there have been some interesting 
attempts to study radiative corrections within its framework [3-5]. To our 
knowledge however, there has never been a complete calculation performed. 
We start our systematic one loop level study of the model from a low-energy 
muon decay calculation. The subject has already been explored qualitatively 
in [6]. The main result of this paper was to show that the quadratic top 
mass dependence of the obhque corrections to Ar is lost. In the SM these 
corrections come from constraints imposed by the SSB sector on the Weinberg 
angle counter-term. Here similar constraints connect the counter-term with 
the heavy SSB scale. As a result the top quark is effectively massless. By 
the same the SM one loop corrections do not constitute a subset of the full 
contributions. Therefore in general, it is not true that one can properly fit 
New Physics Models (NPM) by taking one loop SM corrections modified with 
tree level NPM couphngs. These issues are further explored in [7] . 

The question which we wish to answer is the following: can we or can we not 
accommodate the present experimental life-time of the muon within a model 
that has a minimal Higgs sector structure supported by phenomenology and 
the smallest possible number of unknown free parameters? It is common wis- 
dom that when there are many free parameters any data can be fitted. We 
show however, that this induces a strong correlation between the heavy pa- 
rameters. In fact a full decoupling is not observed, and if the additional masses 
tend to infinity independently, a huge correction results, which is incompatible 
with data. 

We first discuss assumptions on model's structure and parameters. The renor- 
malization scheme is then introduced and the corrections to muon decay enu- 
merated. Numerical estimates follow with a study of the dependence on the 
heavy masses and the heavy symmetry breaking scale. Conclusions together 
with an outlook close the paper. An appendix gathers our notational conven- 
tions and main components of the model. 



2 Structure and Pcirameters of the Model 

As noticed in the Introduction, there are many versions of LRSM. A complete 
analysis at the one loop level requires the model to be fixed. We choose to 
explore the most popular version of the model with a Higgs representation 
with a bidoublet $ and two (left and right) triplets Ai^fi [8]. We also assume 
that the VEV of the left-handed triplet A^, vanishes, (A^) =0 and the CP 
symmetry can be violated by complex phases in the quark and lepton mixing 
matrices. Left and right gauge couplings are chosen to be equal, ql — Qr- 
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We call this model, the Minimal Left-Right Symmetric Model (MLRSM). The 
necessary definitions can be found in the Appendix (for details, see [6,8-10]). 



We also take advantage of several approximations which come from phe- 
nomenological studies. 

(1) Mixing of Fermions 

As usual in one loop analyses, we neglect quark mixings. The case 
of neutrino mixings deserves however additional comments. The effec- 
tive light neutrino mass matrix M — —MdM^^MJ) (Eq. 41) yields three 
light Majorana neutrinos which are predominantly composed of the usual 
"active" neutrinos with a very small 0(l/mjv) admixture of "sterile" 
neutrinos i/£ = vr. The diagonahzation of Mr (Eq. 43) produces 3 heavy 
Majorana neutrinos which are mainly composed of vr. In order to get 
light neutrino masses at the eV scale (as concluded from experimental 
data [11,12]) and without extra relations between Md and Mr matrices, 
niN must be large tun > 10^(10^^) GeV for the lepton (quark) see-saw 
mechanism. However, we would also like to explore lower scales with vr 
of the order of TeV. A crude hght-heavy (LH) mixing estimate 0(l/mAr) 
would give in this case larger couplings. However, they would lead to a 
problem with obtaining the light neutrino spectrum, namely, from Eq. 41 
their masses would be much above the eV scale. A fine-tuning of Mj^ 
and/or Mr parameters or additional discrete symmetries must be ap- 
plied for the full neutrino mass matrix to get the proper light neutrino 
spectrum. Therefore it has been argued in [13] that it is not natural to 
obtain large LH mixings for heavy neutrinos at the TeV scale. In accor- 
dance with these arguments we assume here that the light and heavy 
neutrino sectors are disconnected (negligible mixings). In this way, Wi 
couples only to light neutrinos, while W2 couples to the heavy ones. 
and Z2 turn out to couple to both of them [9,10]. This generates automat- 
ically an extended flavour symmetry, where transitions outside of a family 
composed from a lepton, a light and a heavy neutrino are forbidden. 

As lower limit on the heavy neutrino masses we use the direct experi- 
mental limit from the lack of Z — > vN decay, which is mjv > Mz- 

(2) Mixing of Charged Gauge Bosons 

In principle the model allows for mixing of charged gauge bosons. How- 
ever, experimental data analyses give the following conservative upper 
bound on the mixing angle [14,15] 



The tree level contribution to Ar coming from the mixing is proportional 
to 



^1 < 0.013 rad. 



(1) 




Ml 



(2) 



Ml 
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Even if the second charged gauge boson had a mass of the order of the 
SM W, this number would be neghgible compared to the experimental 
value which is of the order of 3%. 

We therefore put C = 0, which also means that = 0. There are several 
advantages of this approximation. First there is no need to renormalize 
the mixing of the gauge bosons. It turns also out, that together with 
the previous approximation on lepton sector mixings, there is no need to 
renormalize the a priori possible mixing between the light and the heavy 
neutrino within a family. At last, the QED contributions to the process 
form a self-contained class as in the SM (see section 4) . 

This model has the nice feature that all the constraints on the right 
handed sector come uniquely from one loop corrections. The tree level 
Wi exchange diagram is not sensitive to the additional gauge structure 
anymore [16]. 
(3) Yukawa Couplings to Charged Higgs Scalars 

The approximations from the preceding two points, leave still a pos- 
sibility for muon decay through one of the charged Higgs scalars Hi'. It 
turns out, that the experimental data on polarized muon decay asym- 
metries are compatible even with a decay through scalar currents only. 
However, inverse muon decay bounds these contributions to be at most 
one order of magnitude smaller than the SM left-handed current decay. 
We assume in this work that these diagrams are either negligible or re- 
quire only to be included at the tree level in which case, the space left 
for Ar in Eq. 18 would be respectively smaller. 



3 Renormalization in One Loop Order 

As the basic set of input parameters we choose the electromagnetic coupling 
constant and the masses of the four gauge bosons, Higgs scalars and fermions. 
It turns out that as long as no corrections need to be included to tree level 
Higgs scalar exchange diagrams, the on-shell conditions of gauge bosons only 
suffice to fix all of the necessary counter-terms. Moreover, decoupling effects 
should be automatically included. 

In the present approximation, where we neglect charged gauge boson mixing, 
only the Weinberg angle requires renormalization. We recall here its counter- 
term [6] 



Sslr = 2ci 



, {SMl+SMl)-{6M^^+6M^ 
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2 ((Mi + M|J-(M^^+M^j)' 

1 {2Ml + M|J^M|^ + (2M|^ + Ml)6Ml 

2 ((M|,+M|J-(M^^ + M^J)' 



As discussed in section 2, no fermion mixing renormalization is needed, and 
the hard corrections (factorized weak contributions) are properly included 
through the simple fermion counter-terms 

SZ^f!:n = ^'^L,R, (4) 

where we used the following decomposition of fermion self-energy 

S = pP^S^i + pPrJ:^r + Pl^l + Pr^r. (5) 



An interesting problem is connected to charge universality and renormalization 
of the electromagnetic couphng. We wish here to show that charge universality 
follows simply from Ward identities and the constructive proof gives also the 
correct counter-term. To one loop the potentially problematic contributions 
come from diagrams involving a heavy neutrino and the traditional approach 
gives a result independent of these masses only after summation of vertex and 
external line contributions [6]. 

Let us start from the following relation, which comes from charge assignments 
within a fermion doublet and the definition of the physical fields 







amp. 
^2 I amp. 
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amp. 



, (6) 



where C/q is the bare neutral sector mixing matrix Eq. 31 multiplied by the 
renormalization constants of the physical fields 
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and (. . .)amp. is a shorthand for amputated Green functions. Prom this we 
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obtain 



After taking the divergence of the current, we can use the U{1) Ward iden- 
tity for the B field and the on-shell renormahzation conditions on the fermion 
propagators and the electromagnetic vertex, which leads to the following iden- 
tity 

e = ; '° „ (^0^)33, (9) 
^cos2e?i. 



which can be put into the following form 



1 



en 1 sin 0° tan — cos 6^ J cos 2Q9^ tan ^ 
^(Z^ + "J 5^ 5^Z|„ ( 10) 

e i/cos2eyi. 



cos0°tane^ -Fsin0°Jcos2e^tane^ _i 



cos 26^ 



Since none of the above renormahzation constants depends on the initial 
fermion species, we have obtained the required charge universality. At the 
same time we can expand this relation to first order to yield the electromag- 
netic coupling renormahzation counter-term 



5e 1 

(11) 

sin (f) tan — cos 0-\/cos2©ty tan ^\ 
Vcos2eH. 

cos tan ©vt^ -I- sin 0-\/cos2©^ tan ©^^ i 
^ Vcos2©H. 

We have checked by explicit calculation that the above formula gives the same 
value as the usual approach. We would like to stress that to our knowledge, 
such a formula for LR models has never been derived, although similar meth- 
ods have been used in SM analyses [17]. 



4 Structure of Corrections to Muon Decay 



The muon life-time is parametrized through the Fermi coupling constant, the 
mass of the muon and the QED corrections to the four-fermion interaction Aq, 



6 



which are presently known up to second order in the fine structure constant 
[18] 



1 _ Glml 
r„ 1927r3 



(1 + A?). 



(12) 



The Fermi constant on the other hand is related to the tree level SM coupling 
of the charged W boson to fermions through 



where Ar are higher order corrections to which we already made reference. 

There are two problems with these formulas, when moving from the SM to 
other interactions. Let us first consider Eq. 12. It is based on the assumption, 
that the basic process is described by the four-fermion interaction, which in 
the charge conserving form is of a pure V — A type. In fact, as long as the 
interaction has only an admixture of vector and axial currents, the QED 
corrections are finite and gauge invariant, hence meaningful. Notice however, 
that if the process is induced also by right handed currents, then after moving 
to the charge conserving form of the interaction (Fierz transformation), there 
appear also scalar and tensor interactions, which are known not to have a 
finite QED correction. The same problem occurs if we add charged scalar 
particles to the list. There are two possibilities to remedy the situation, the 
first being of course calculating by any means the process in the full model 
and resign from the separation of QED corrections. The second possibihty 
is somewhat simpler. If the tree level corrections from right-handed and/or 
scalar interactions are of the same size as the one loop corrections to the basic 
diagrams, we can simply ignore one loop contributions to these additional 
currents and consider Eq. 12 as approximate and valid to one loop order only. 

The second problem we have to face is the fact that the tree level coupling to 
the light charged gauge boson can be different from the SM one. This concerns 
mainly the sine of the Weinberg angle sw. In fact this happens to be the case 
of the considered model, where due to constraints if we fix the mass of the two 
light gauge bosons, then sw is given by a function of the heavy SSB scale vr. 
This dependence is depicted in Fig. 1. For small values of vr, the difference 
from the SM value is large. We choose here to include the change of sw from 
the SM to the LRSM in Ar. 

Ar can now be obtained from the formula 



Gf 



(1 + Ar), 



(13) 



^/2 SM^sl 
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Ar= (11^(^5^21^ + 2^-^ + .. + ..) (14) 

{Sw)LRSM V e stv J 

_ {s\y)lrsm — {S]y)sM 
{sw)lrsm 

where Sy denotes the vertex corrections, which consist of the proper one loop 
vertex diagrams and the incomplete counter-term made only of the fermion 
wave function renormalization constants 

Sv = ^^{K..w + ^,.,w) + \{5Zl + SZ'i: + 5Z^ + 5Zl^), (15) 

with A being the coefficient in front of the operator ^'^Pl, and 5b represents the 
box contributions. The last term comes of course from the "renormalization" 
of the Weinberg angle between the two models with {sy^)sM = 1 — -j;^ and 
{^w) LRSM as obtained by solving Eqs. 28,29. 

The strong dependence on the light fermion masses in 5e is avoided as usual 
by a shift up to the Zi mass, and insertion of the running of the fine sturcture 
constant, for which we take [20] 

Aa{Mz,) = 0.059394 ± 0.000395 (16) 

The factorization of the QED corrections is obtained with the Sirlin's method 
[19], which amounts to rejecting the infrared divergent box diagram and re- 
placing the photon vanishing mass by the Wi mass in the infrared divergent 
lepton wave function renormalization constants. 



5 Quantitative Results 

The evaluation of one loop corrections within the LRSM is a task of moderate 
size as far as the number of diagrams is concerned. In fact approximately 600 
had to be calculated already after our simplifying assumptions. It would not 
be possible to perform this work without using an automated system. For the 
generation of diagrams we used the C-|— I- hbrary DiaGen [21] , which currently 
contains a topology generator, with several tools to analyse the properties 
of the created objects, and a diagram generator with support for Majorana 
fields. The output has then been algebraically simplified with FORM [22], 
and at last numerically evaluated with the help of the FF [23] based package 
LoopTools [24]. 
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As discussed in the previous section we parametrized the muon hfetime cor- 
rections coming from the LRSM through Ar, which is defined analogously as 
in the SM. With the present values of the coupling constants and masses [15] 



Gf = 1.16639(1) ■ 10"^ GeV-^, l/a. = 137.0359976 ± 0.00000050, 
Mw, = 80.451 ± 0.033 GeV, Mz, = 91.1875 ± 0.0021 GeV, 



(17) 



the value of Ar with error is 



Ar = 0.032 ±0.004. 



(18) 



We depicted this experimentally allowed range by a shaded region on the 
relevant figures. 

As noted already in a previous work [6], we should not expect decoupling in 
the sense that for large Vji and large masses of the additional particles the 
SM result for Ar would be obtained. Some type of decoupling is however 
observable. For example if we take the box diagrams in the 't Hooft-Feynman 
gauge, then the result tends to the SM one as depicted in Fig. 2. It is worth 
noting however the effect of taking heavy neutrinos with a low vr as for the 
(c) curve, where the contribution blows up. This is simply a consequence 
of the fact, that the ratio of a neutrino mass and the heavy SSB scale is 
proportional to the Yukawa coupling Hm (Eq. 36) and the respective diagrams 
are proportional to at least the square of these couplings. Obviously, if the 
Yukawa couplings start to be larger than one then the perturbative expansion 
must break down. 

An interesting effect is obtained, if we take the masses of the Higgses to follow 
some simple pattern as in the Appendix Eqs. 22 and 23. The respective Ar is 
shown for several heavy neutrino masses in Fig. 3. With growing vr the value 
grows strongly away from the allowed range and these parameters must be 
rejected. Although heavy neutrinos lower down Ar, we cannot obtain a rea- 
sonable value even if their masses are at the edge of the perturbatively range. 
The line (d) realizes this situation with the largest possible heavy neutrino 
mass as a function of vr (Eq. 36) 



If we now assume for simplicity that all of the Higgs scalar masses are equal, 
apart from the SM Higgs boson, then we obtain the strong dependence as 
depicted in Fig. 4. If all the scalars arc approximately two times heavier than 
(for large Higgs masses), the experimental value for the muon decay life- 
time can be accommodated. Let us note at this point that large Higgs masses, 
at least of the order of a few TeV are needed because of FCNC [4] . It is obvious 



(19) 
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from Fig. 4 that Higgs scalars, heavy neutrinos and additional gauge boson 
masses are very much fine-tuned to be within the SM gray area, e.g. the hne 
(a') with ruH = 1 TeV and tun = 100 GeV gives vr ~ 800 GeV, which fixes 
and M^j to the values as in Fig. 5. 

It is interesting, that for a larger SSB breaking scale, the variation of Ar with 
the SM Higgs scalar mass is negligible. This is shown in Fig. 6 for vr = 2390 
GeV and neutrino and heavy scalar masses chosen to fit the experimental 
value. 

At last let us comment on the dependence on the top quark mass. We show 
the contributions of the third quark family in Fig. 7. The values of the top 
mass spread over a vary large range to show the behaviour of the correction. 
As already forseen in [6], for low vahics of Vr the variation is described by a 
negative quadratic function. However for a Vr as low as 1 TeV, only a positive 
logarithmic contribution is visible. Notice also, that even in the low vr range, 
the top mass squared enters with a smaller coefficient than in the SM. 



6 Conclusions 



In this paper we have studied the full one loop corrections to the muon decay 
in a self consistent Left-Right symmetric model. We have shown quantitatively 
that the contributions have a different structm^c from the SM ones and that 
they cannot be separated into these and some corrections that would vanish 
with Vr. Moreover, we have shown that the muon decay alone already puts 
some stringent restrictions on the different heavy particle masses with respect 
to the heavy SSB breaking scale. 

Our analysis should be extended to cover also other low-energy experiments 
[14]. This should elucidate the question of the contribution of the boson 
to the muon decay. It will then also be possible to derive bounds on the extra 
boson masses. 

At last let us note that several of the assumptions that we here took could be 
raised, but it is doubtful that this would change qualitatively the numerical 
results. On the other hand it would certainly make the analysis more involved, 
starting from the necessity to renormalize the ^ angle, LH neutrino mixing and 
ending with problems with the QED contributions. 
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7 Appendix 

In this appendix we gather our definitions and give a short account of the 
particle content of the model. 

7.1 Higgs sector 

The Higgs sector contains one bidoublet and two triplets 



of which K2 and vl are assumed to vanish. 

The full potential has been studied in [8-10]. Here we only recall the physical 
spectrum of the particles, which consists of 

(i) four neutral scalars with J^^ = 0++ i = 0, 1, 2, 3) , 

(ii) two neutral pseudoscalars with J-^^ = 0+" {A^ i — 1,2) , 



If vr ^ Ki and all the parameters of the potential which enter the Higgs 
masses are taken to be 1 (these are combinations of /ii,2, Ai,...,6, pi,...,4 defined 
in [8]), then neglecting terms proportional to the VEV of the SM, the masses 
satisfy the relations 



Mh^ = Mjjo = Mho = Mao = M^o = = = = vr/V2, 




(20) 



with the allowed Vacuum Expectation Values (VEV) 




(21) 
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(22) 
(23) 

(24) 



7.2 Gauge boson sector 



Gauge boson masses are generated by the following mass terms (with the 
assumption of equal couplings for the two SU{2) groups) 









+ h. 










V 





with 



^^^Charged ' 



9_ 

4 



+ 



-2KiK2 



-2kiK2 k+ + 2vr ' ' 



and 



(25) 



(26) 



Neutral 



2 



-44 4(4 + 44) 



V 



-2gg V 



'n.2 

R 



-2gg'v^ 
2g% 



R 



R 



(27) 



where k,+ — \Jk{-\- kI- The masses of the physical gauge bosons are then given 

by 



- — 



^i.. 4 



(28) 
(29) 
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The symmetric mass matrices are diagonalized by the orthogonal transforma- 
tions 



V 



-sin^ cos^ 



(30) 



and 



V ^ J 



( 



CwC Cw-S Sw 

—swSmc — Cms —swSms + cmc cwSm 

— SwCmC + SmS —SwCmS — SmC CwCm 



(31) 



where 
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sinQw ' 



9 = 



■\/cos26vK 



, cw — cos Qw, — sin Qw, 



Vcos20^ 

cm = 7^ , SM — tgK^w, s — smcp, c — coscp. 

cos <dw 

The mixing angles are given by 



tan 2^ = 2 — ) sm20 = 



g'^Kl_^/ cos 2Qw 



2cos'Qw {mi -Ml)' 



(32) 



Fig. 5 sums up the mass dependence of the additional gauge bosons and two 
sets of Higgs scalar particles Eqs. 22,23 on the vr scale. 



7.3 Neutrinos 



The MLRSM naturally contains left as well as right handed neutrino states. 
We chose the following basis for these fields 









(.A 










riR = 

















(33) 



where both uk^r) form 6-dimensional vectors. The neutrino mass matrix takes 
the form 
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^ Md^ 



Mr J 

Mn = -^{hiKi + hiK2) = M]), 



Mr = V2hMVR = Ml. 

hi, hi, hM are the respective Yukawa coupling matrices. 

Ml, can be diagonalized with the following unitary transformation 

(kt\ 



(34) 

(35) 

(36) 



U 



L 



(37) 



where the K^ r matrices have dimension 6x3. The LEP neutrino counting 
results show that there must be three hght active neutrino states. This means 
that we must have M^ ^ Mr, which after requiring "natural" couplings 
(order one), turns into ki^2 <^ vr. 

Let us now introduce 3x3 matrices Ua{h), Urkj^) [9]: 



Kr^ 



(38) 
(39) 



The diagonalization equation assumes the form {nidiag and M^iag correspond 
to light and heavy neutrino mass matrices, respectively): 



u 



LI ^ m 

Lh 



Ml 



Md 

Mr 



TT* TT* 
Uri Urh 



This can be written as: 



/ mdiag 
V Md^ag 



(40) 



Ul (-MnM^'Ml)uii^maiag, 

Um^-M^'MlUli, 

Ul^MRURh^Mdiag, 

ui,c^mum*r)-'Ur,, 



(41) 
(42) 
(43) 
(44) 
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where the unitarity of U, and the large scale difference Mdiag ^ rndiag, have 
been used. Two important conclusions can be drawn from it: 



• the matrices Uli and Unh are approximately unitary, 

• the elements of the non-diagonal submatrices Uri and Uih are small, of 
order ^ < ^^(i^, where = {M^^ag). 

The symbol (...) denotes the relevant scale of mass matrices. 

Altogether we can write 

^-f5)-f.,?<i*. °'±"*V (45) 



VA-;,; vo(i/m«) 0(1) 
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Fig. 1. sin^ @\Y as function of vr with M\Yi and Mz-^ as in Eq. 17. 
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Fig. 2. Contribution of box diagrams to Ar. With increasing vr heavy particles 
decouple and the lines aim at the SM contribution. The (a) line is for (three heavy 
neutrinos) mjv = 100 GeV; (b) is for rriM = 500 GeV; (c) is for m^v = 2 TeV. Higgs 
particle masses obey Eqs. 22,23. The gray area shows the experimentally allowed 
values of Ar. 
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Fig. 3. Ar as function of vr for different heavy neutrino masses. Higgs masses are 
chosen according to Eqs. 22,23. The (a) hne is for (three heavy neutrinos) mj^ = 100 
GeV; (b) is for mjv = 500 GeV; (c) is for mjv = 2 TeV. Line (d) shows the results 
when heavy neutrino masses follow from Hm = 1 (see Eq. 19). The gray area shows 
the experimentally allowed values of Ar. 



18 



0.2 




500 1500 2500 3500 

[GeV] 



4500 



Fig. 4. Ar as function of vji. Sets with and without primes show results for three 
heavy neutrino masses with ttin = 100 GeV and ttin = 2 TeV respectively. The 
lines describe different values of Higgs scalar masses: (a) is for all Higgs masses 
Mh = 1 TeV; (b) is for Mh = ^ TeV; (c) is for Mh = 10 TeV. The gray area shows 
the experimentally allowed values of Ar. 




Fig. 5. Masses of additional gauge bosons and two sets of Higgs scalar particles 
Eqs. 22,23 as function of the vr scale. 
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Fig. 6. Ar as function of the lightest Higgs scalar mass M^o.The gray area shows 
the experimentally allowed values of Ar and the heavy particle spectrum is chosen 
to fit approximately to this region, namely, vr = 2390 GeV, rriN = TeV, mn = 5 
TeV. 
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Fig. 7. The contribution of the third quark family to Ar as function of vr for 
different top quark masses. 
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